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Control of synchronization in coupled neural systems by time-delayed feedback 
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We discuss the synchronization of coupled neurons which are modelled as FitzHugh-Nagumo sys- 
tems. As smallest entity in a larger network, we focus on two diffusively coupled subsystems, which 
can be interpreted as two mutually interacting neural populations. Each system is prepared in the 
excitable regime and subject to independent random fluctuations. In order to modify their coop- 
erative dynamics, we apply a local external stimulus in form of an extended time-delayed feedback 
loop that involves multiple delays weighted by a memory parameter and investigate if local control 
applied to a subsystem can allow one to steer the global cooperative dynamics. Depending on the 
choice of this new control parameter, we investigate different measures to quantify the influence on 
synchronization: ratio of interspike intervals, power spectrum, interspike interval distribution, and 
phase synchronization intervals. We show that the control method is more robust for increasing 
memory parameter. 



1. INTRODUCTION 

The network of neurons in the brain exhibits a sub- 
tle balance of dynamic chaos and selforganized order. A 
number of neurological diseases like Parkinson or epilepsy 
are characterized by a disturbance of this balance, e.g. 
synchronized firing of electrical pulses of the neurons 
Schiff et"all . Il994{ Modern concepts of time-delayed 



feedback control have recently been appl i ed to sup- 



press this undes ir ed syn c hrony iGassel et a/.l. I2007ll200g : 
PoDOvvch eTai\. l2005l l2006l: iRosenblum k Pikovskv . 



2004aU bl: ISch511 fc Schuster!, l200, 

It was shown earlier jHauschildt et 'd\, l2006j that 
application of time-delayed feedback, which was origi- 
nally suggested for deterministic chaos control Pvragasl 
1992{ . is able to influence the cooperative dynamics. 
As a measure of cooperative dynamics, we consider 
coherence, timescales, and synchronization of noise- 
induced oscillations. With this method, a control force 
is constructed from the difference of the current state 
of a system to its time-delayed counterpart. Previ- 
ously, it has also been used to influenc e noise-induced 
oscill a tions of a sing l e exci t able system iBalanov et al 



2004: Janson et al 



Prager et all l2007l 



tion 



|2004 iPototskv fc Jansonl . 12005 : 
of s ystems below a Hopf b i furca- 



Flunkert fc Schdlll. l2007t iPomplun et all l2005l 



iPototskv fc Ja nson'. '20 071: ISchoU et a /. 2005*1 or a globa l 



bifur cation 
l2008t . and of 



Hizanidis et al. 



systems 



spatially 
Balanov et all 



IStegemann et al. . 2006l | 



20061 .Hizanidis fc Scholj 
exten de d reaction- dif fusion 
l2006l : ISchoU et all l2008l 



In t his paper, we extend the work of iHauschildt et oil 
2006j | by application of a different feedback stimu- 
lation with multiple time delays. External delayed 
feedback loops have been sug gested for supp r ession 
of pathological brain rhythms Popovvch et all l2006t 



* Elec tronic address: |schoell@physik.tu-berIin.de| 

URL: I http : //www ■ Itp ■ tu-berlin . de/schoell | 



iRosenblum fc Pikovskvl l2004bt . Our method, also 
known as extended time-delayed feedback, was initially 
proposed by Socolar et al. in order to extend the 
domain of effective s tabili zation of unstable periodic 
orbits [Socolar et all Il994| . It generalizes the Pyra- 
gas scheme by introducing an additional memory pa- 
rameter and is known for successful stabilization at a 
larger r ange of parameters compared to the Pyragas 



larger r ange or pararneters compared to tne ryragas 
method Dahms et al 1. 120071 120081 : f Schlesner et al. 1. 120031 



lUnkelbach aZ.I . l2003l | . In case of noise-induced oscilla- 
tions, extended feedback has been demonstrated to result 
in drastically im proved coherence and arbitrarily long 
correlation times [Pomplun et aLl . [2007 |. 



Our aim is to control the global cooperative dynamics 
of the ensemble of neural populations by local application 
of a stimulus to a single subsystem that involves multiple 
delays in the feedback. In particular, we are interested in 
the study of the effects of the memory parameter as new 
control parameter on the synchrony properties in coupled 
neural systems. We compare this to the case of vanishing: 
memo ry parameter, as investigated by IHauschildt et al\ 
200611. and flnd that the control scheme is more robust 



for multiple time-delayed feedback in the sense that en- 
hanced synchronization is independent on the tuning of 
the time delay. 



The structure of the paper is the following: In Sec. [51 
we introduce the model equation and the control method. 
In Sec.[3l we discuss the configuration of the uncontrolled 
system and consider the average interspike intervals. Sec- 
tions m and [S] are devoted to the effects of the control 
parameters on the power spectrum and the interspike in- 
terval distribution, respectively. In Sec. [6l we introduce 
a phase variable and discuss effects of phase synchroniza- 
tion. Finally, we conclude with Sec. \7\ 
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2. MODEL 



111 the following, we consider two mutually coupled 
neurons modelled by FitzHugh-Nagumo systems: 

eiii(0 = x,{t) - ^ - y,it) + C[x2it) - x^it)] 
yi{t) = xi{t) + a + D,^,{t) (la) 

OO 

+K J2 [yi {t-{n + l)r) -yi{t~ nr)] 



x^(t] 

e2X2{t) = X2{t) - - y2it) + C[xi{t) ~ X2it)] 



mit) = X2{t)+a + D2i2{t), 



(lb) 



where Xi, yi and X2, y2 correspond to single excitable sys- 
tems representing neurons or neural populations, which 
arc diffusively coupled in the activator variables xi, X2 
with coupling strength C. The variables yi and 2/2 rep- 
resent the inhibitor. Throughout this paper, we consider 
neurons in the excitable regime at which no autonomous 
oscillations occur. Thus, we fix the excitability parame- 
ter a as a = 1.05. 

In order to introduce different timescalcs in the two 
subsystems, we choose ei = 0.005 and 62 = 0.1. Each 
neuron is driven by Gaussian white noise S,i{t){i = 1,2) 
with zero mean and unity variance. The noise intensi- 
ties are denoted by parameters Di and D2, respectively, 
where we keep D2 fixed at D2 = 0.09 in the following. 
The last term in Eg. ITal) describes ext ended time-delayed 
feedback control Socolar et all Il994l | with time delay r, 
feedback gain K, and memory parameter R. Note that 
there exists an equivalent recursive form of the feedback 

F{ty. 

00 

F{t) = Kj2R''[yi{t-in + l)r)-yi{t-nT)] (2) 



?i=0 



K[yi{t-T)-y,{t)]+RF{t-T). 



(3) 



The latter form is more amenable to experimental real- 
ization and for practical applications because the delayed 
feedback force F{t — r) replaces the infinite series. 



2.1. Moderate, weak, and strong synchronization 

The two neurons or neural populations of Eqs. ([T]) are 
prepared in the excitable regime. Without external in- 
put, they remain in their stable fixed points. Random 
fiuctuations lead to spiking. 

For reason of comparison, we consider at first the case 
where no control is applied to the system. We set the 
noise intensity D2 in the second subsystem to a small 
value, D2 = 0.09, to reafize some background noise level. 
Depending on the coupling strength C and the noise in- 
tensity Z^i, the two subunits show cooperative dynamics. 

Calculating the average interspike interval of the two 
neurons (Ti) and (T2) and their ratio, one can see 
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FIG. 1: (Color online) Ratio of interspike intervals {T\) / {T2) 
of the two subsystems in dependence on the coupling strength 
C and noise intensity D\. No control is applied to the sys- 
tem. The dots mark the parameter choice for different syn- 
chronization regimes used in the following. Other parameters 
as in Fig. [2] 



in Fig. [T] how the frequency synchronization changes 
in dependence on the coupling strength and noise in- 
tensity. For a small value of Di and large coupling 
strength, the two subsystems display synchronized be- 
havior, (Ti) / {T2) « 1. On average, they show the same 
number of spikes indicated by bright (yellow) color. 

Since we are interested in the effects of a control force 
on the synchronization, in the following, we consider 
three different cases: moderately, weakly, and strongly 
synchronized systems given by the specific choices of the 
coupling strength and noise intensity in the first subsys- 
tem (C = 0.2, Di = 0.6), (C = 0.1, Di = 0.6), and 
(C ~ 0.2, Di = 0.15), respectively. These different cases 
of stochastic synchronization are marked as black dots in 
Fig. [TJ Strong synchronization can be found for small 
noise intensity Di and large C, e.g., Di ~ 0.15 and 
C ~ 0.2. Moderate synchronization is given for a choice 
of Di = 0.6 and C = 0.2, and weak synchronization can 
be realized by Di = 0.6 and C = 0.1. 

In terms of the time series, the three configurations 
of moderate, weak, and strong synchronization are dis- 
played in Fig. [2] as panels (a), (b), and (c), respectively. 
In all panels, the black, grey (red), and lightgrey (green) 
curves refer to the global xs-variable, the X2-, and xi- 
variable, respectively, where xs is given by xy: = xi+X2- 

Moderately synchronized systems perform mostly syn- 
chronized spiking. However, there are certain events 
where only one system shows an oscillation (see panel 

(a) of Fig. [21). In the case of weak synchronization, the 
spikes of the two subsystems coincide less as can be seen 
from the time series of the summarized signal in panel 

(b) of Fig. [21 For strongly synchronized subsystems (see 
panel (c) of Fig. the time series of the xi- and 2:2- 
variable exhibit spiking at the same time. 
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FIG. 2; (Color online) Time series of the coupled FitzHugh- 
Nagumo system in the absence of control. Panels (a), (b), and 
(c) correspond to moderately (C = 0.2, Di — 0.6), weakly 
(C = 0.1, Di = 0.6), and strongly (C = 0.2, Di = 0.15) syn- 
chronized systems, respectively. In all panels, the black, grey 
(red), and lightgrey (green) curves refer to the summarized 
variable xs = xi + X2, the X2-, and the Xi-variable, respec- 
tively. Other parameters: a = 1.05, ei = 0.005, £2 = 0.1, and 
D2 = 0.09. 



2.2. Time series and control 



Equation (|lap already includes the time-delayed feed- 
back scheme, whose parameters are the feedback gain 
K, the time delay t, and the memory parameter R. 
The spe cial case of R = 0, also known as Pyragas 
control iPvragas , Il992ll , was already investigated by 
2006[. Therefore, in the present work, 



Hauschildt et al. 



we discuss an extension to a different feedback stimula- 
tion involving multiple delays. In this extended form, the 
method generates a control force from the differences of 
the states of the system that are one time unit r apart. 
The memory parameter R G (—1,1) can be understood 
as a weight of states that are further in the past. For 
vanishing memory parameter i? = 0, only one delayed 
state enters the generation of the control force. 

We consider the case where the control force is ap- 
plied to the first system only and in the inhibitor vari- 
able. Note that one could also realize a feedback scheme 
that is applied to both subsystems and study effects of 
different values of the control parameters for each sub- 
system, but these investigations are out of the scope of 
this work. Here we put special emphasis on the effects 
due to changes in the memory parameter. 

Figure [3] depicts different time series for moderate syn- 
chronization (C= 0.2, Di — 0.6) in the absence of con- 



FIG. 3; (Color online) Time series of the coupled FitzHugh- 
Nagumo system for moderate synchronization. Panel (a) cor- 
responds to no control. In panels (b) and (c), time-delayed 
feedback is applied to the system with different memory pa- 
rameters i? = and R = 0.9, respectively. Other control 
parameters are fixed at r = 1 and K = 1.5. Grayscale (color 
coding) and other parameters as in Fig. (2] 



trol (panel(a)) compared to the cases when time-delayed 
is applied to the system (panels (b) and (c) for memory 
parameters of i? = and R = 0.9). The time delay and 
feedback gain are chosen as r = 1 and if = 1.5. In all 
panels, the black, grey (red), and lightgrey (green) curves 
correspond to the summarized global signal xs, the X2-, 
and the xi-variable, respectively. One can see that the 
time-delayed feedback enhances the synchronization, i.e., 
events at which only one system oscillates are less fre- 
quent. In this sense the choice of i? = is more efficient 
compared to larger memory parameters. 



3. RATIO OF AVERAGE INTERSPIKE 
INTERVALS 

As first measure to quantify changes in the synchro- 
nization due to the control force, in this section, we 
consider the ratio of average interspike intervals. In 
the presence of a control force, i.e., if 7^ 0, the coop- 
erativity can be influenced by the control parameters. 
Varying the feedback gain if, the time delay r, and 
the memory parameter i?, the average interspike inter- 
val c an be altered. The case of i? = was discussed in 
Ref. [Hauschildt et aLl . l2006t . 



For fixed feedback gain if 1.5 and two different val- 
ues of i? (i? = and R = 0.9), Fig. [J depicts the av- 
erage interspike intervals of the subsystems, shown as 
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FIG. 4: (Color online) Interspike intervals in dependence on the time delay r. The green dots correspond to the ratio of 
the interspike intervals {Ti) / {T2) of the two subsystems, which are also depicted by solid (black) and dashed (red) curves for 
(Ti) and {T2), respectively. Panels (a), (b), and (c) correspond to the case of moderate, weak, and strong synchronization, 
respectively. The control parameters are chosen as, K = 1.5 and the memory parameter corresponds to 7? = and 0.9 in the 
top and bottom panels, respectively. Other parameters as in Fig. [21 



solid (black) and dashed (red) curves for (Ti) and (r2), 
and their ratio (as dotted (green) curve) for the case of 
moderately, weakly, and strongly synchronized systems, 
respectively, in dependence on the time delay r. The left, 
middle, and right panels correspond to the case of mod- 
erate, weak, and strong synchronization, respectively. 

In all cases, tuning of the time delay leads to either 
enhanced and deteriorated synchronization for i? = 0. If 
more information of past states (i? ^ 0) is included, how- 
ever, the variation of the ratio (Ti) / {T2) is less sensitive 
to the specific choice of r. The bottom panels of Fig. 3] 
do not show large deviations for the ratio of average in- 
terspike intervals. Therefore, a larger memory parameter 
renders the control method more robust. 

One can also consider two-dimensional projections of 
the control parameter space spanned by if, t, and R. 
Parameterized by the feedback gain and time delay, Fig. [5] 
displays the ratio of (Ti) and (T2) for moderate, weak, 
strong synchronization, respectively. The four panels in 
each figure correspond to a memory parameter of i? = 0, 
0.35, 0.7, and 0.9. Note that Fig. [4] can be understood as 
horizontal cuts for K = 1.5 in the respective diagram of 
Fig.li 

As also shown in the one-dimensional projections (see 
Fig- HI) J £^1^ increase of the memory parameter R at fixed 
K yields smaller changes of the ratio (Ti)/(T2). Indepen- 
dent on the feedback gain K, the desynchronized darker 
region in Fig. [5] at a time delay t « 3 is much less pro- 
nounced. For R = 0.9, the ratio of average interspike 
intervals is constant in a wide range of the whole control 
domain, which reflects the robustness of this extended 
feedback method. 



4. POWER SPECTRUM 

The two uncontrolled neural populations have differ- 
ent timcscalcs, i.e., ei 7^ €2, and therefore, oscillate at 
different frequencies. We have shown in the last section 



that the mean frequencies of the subsystems measured 
by the average interspike intervals match for certain con- 
trol parameters. In order to gain a better understanding 
of the shift in the timcscalcs in the subsystem in the 
presence of control, we investigate the power spectrum 
in this section, where we focus on the role of the memory 
parameter. 

At first, we consider the case when no control is ap- 
plied to the system, i.e., K = 0. The result can be 
seen in Fig. [6] where the solid (red) and dashed (blue) 
curves correspond to the power spectrum, of the xi- and 
a;2-variable, respectively. The left, middle, and right 
panel refer to the case of moderate, weak, and strong 
synchronization, respectively. This change of frequency 
synchronization is reached by choosing different coupling 
strengths C and noise intensities Di. In the case of mod- 
erate synchronization (C = 0.2, Di = 0.6) the peaks 
of the power spectrum partly overlap. Weak synchro- 
nization (C = 0.1, Di = 0.6) shows a smaller overlap, 
whereas for strong synchronization (C = 0.2, Di = 0.15) 
the spectra almost coincide. 

Next, we apply extended time-delayed feedback to the 
system. Figures [3 [51 and M show the power spectrum 
of the summarized signal xs = xi + X2 for the three 
above mentioned cases of synchronization in dependence 
on the time delay r. In all figures, panels (a), (b), and (c) 
correspond to a memory parameter of i? = 0, i? = 0.5, 
and R ~ 0.9, respectively. The feedback gain is fixed at 
K ^ 1.5. 

It can be seen that, depending on the choice of r, the 
main frequency component is shifted. Thus, the con- 
trol scheme is able to support different timcscalcs. Note 
that, for instance in the case of moderate synchroniza- 
tion (Fig. I?}, the control force enhances the frequency 
corresponding to the dynamics of the xi- or a;2-variable. 
Compare the bright (yellow) areas in Fig. [Tito the middle 
panel of the uncontrolled case in Fig. [SI A time delay of 
T w 3 favors a frequency of / ~ 0.4, i.e., the main fre- 
quency of I'l, and r w 5 enhances components of / « 0.2, 
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FIG. 5: (Color online) Ratio of interspike intervals {Ti) / {T2) in dependence on the feedback gain K and the time delay r for 
moderate (C = 0.2, D\ — 0.6), weak (C = 0.1, Di = 0.6), and strong (C = 0.2, Di — 0.15) synchronization in the left, middle, 
and right panel, respectively. The memory parameter are fixed ar i? = 0, 0.35, 0.7, and 0.9 in the four subfigures of all three 
panels. Other parameters as in Fig. [51 




FIG. 6: (Color online) Power spectrum of the two subsystems in the absence of control. The left, middle, and right panels 
correspond to the case of moderate, weak, and strong synchronization, respectively. The solid (red) and dashed (blue) curves 
refer to the xi- and a;2-variable, respectively. Other parameters as in Fig. [J] 



which corresponds to the dynamics of X2- 

For larger memory parameters R, this effect is less pro- 
nounced. The power spectra of the controlled system 
display reduced sensitivity on the specific choice of the 
time delay. The main peak of the spectrum is stronger 
localized at the frequency of the second subsystem, but 
the value of the power spectrum at this main frequency 
is much lower. 

Figure [TO] shows the power spectrum of the global sum- 
marized signal for the three different case of synchro- 
nization and fixed time delay t = 5. In all panels, the 
solid (black), dashed (red), and and dotted (green) curves 
correspond to a memory parameter of R = 0, 0.5, and 
0.9, respectively. The feedback gain is fixed at K — 1.5. 
Panels (a), (b), and (c) of Fig.[TO|can be seen as horizon- 
tal cuts (at r = 5) of Figs. [3 [HI and [HI respectively. 

Opposed to the uncontrolled case (see Fig. [H]), the 
power spectra exhibit a series of distinct peaks that are 
located at the harmonics of the main frequency / « 0.2. 
They are due to the control force which enhances not 
only the main frequency, but also frequency components 
of the higher harmonics. For increasing memory parame- 



ter R, the background becomes broader. We stress that a 
similar effect was found in the context of extended time- 



Pomulun et al. 
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and in a 
where, 



in the Van-der-Pol system 
reaction-diffusion system 
in addition, the peaks become sharper. 



5. INTERSPIKE INTERVAL DISTRIBUTION 

In the previous Section, we have discussed the modu- 
lation of timescales due to time-delayed feedback. To ob- 
tain further information about the timescales present in 
the coupled system, we investigate the probability distri- 
bution of the interspike intervals in the following, where 
we restrict our investigation to the analysis of the sum- 
marized signal = xi + X2- 

Figures [TTJ [T^l and [TO] depict the dependence of the 
distributions of the interspike interval Ts for moderately 
(C = 0.2, Di = 0.6), weakly (C = 0.1, L»i 0.6), and 
strongly synchronized (C = 0.2, Di = 0.15) subsystems 
on the time delay, respectively. These are the three cases 
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FIG. 7: (Color online) Power spectrum of the summarized signal xt, = 3-1+2:2 for moderate synchronization (C = 0.2, D\ = 0.6) 
in dependence on the time delay r. Panels (a), (b), and (c) correspond to a memory parameter oi R = Q, R = 0.5, and R = 0.9, 
respectively. The feedback strength is fixed at K = 1.5. Other parameters as in Fig. [2] 
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FIG. 8: (Color online) Power spectrum of the summarized signal xj: = xi-\-X2 for weak synchronization (C = 0.1, Di = 0.6) in 
dependence on the time delay r. Panels (a), (b), and (c) correspond to a memory parameter of R = 0, R = 0.5, and R = 0.9, 
respectively. Other parameters as in Fig. [7] 



marked in Fig. [T] The corresponding power spectra are 
shown in Figs. [71 [HI and [5] Panels (a), (b), and (c) in each 
figure refer to a memory parameter of i? = 0, 0.5, and 
0.9, respectively. Note that all three figures [TTl [T^ and 
[T51 display the same grcyscale (color code) . 

For increasing memory parameter R, the distribution 
exhibits a weaker dependence on the time delay t. This 
can be seen, for instance, in panels (c) of Figs. [11] and 
[T2I For smaller R values, another effect becomes appar- 
ent. There are a minima in the distribution for inter- 
spike intervals slightly larger than r, indicated by black 
stripes in the {T^,t) plane. These minima arc located 
below the line ~ t (white line in Fig. [TlT a)). A sim- 
ilar, but less pronounced effect can be observed if the 
interspike intervals match integer multiples of the time 
delay as shown by the grey line for Ts = 2r in the same 
panel. We stress that for strongly synchronized subsys- 
tems this structuring of the interspike interval distribu- 
tion becomes more visible as displayed in Fig. 1131 The 
probability distribution becomes multimodal with peaks 



centered near = nr, n = 1, 2, . . . . 

To summarize this Section, the introduction of a large 
memory parameter renders the control method more ro- 
bust against the specific choice of the time delay as is 
depicted in Fig. [T3| for R = 0.99. However, for small 
memory parameters, there is a competing effect which 
structures the distribution in the sense that interspike 
intervals slightly larger than the time delay of the feed- 
back are suppressed. 



6. PHASE SYNCHRONIZATION 

The measures for cooperative dynamics considered so 
far are insensitive to phase relations. A measure of phase 
synchronization can be obtained by introducing phase 
variables for each subsystems, and monitoring their dif- 
ference. 

For definition of such a phase variable, one can gencr- 
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FIG. 9: (Color online) Power spectrum of the summarized signal xt, = x\+X2 for strong synchronization (C = 0.2, Di = 0.15) 
in dependence on the time delay r. Panels (a), (b), and (c) correspond to a memory parameter oi R = G, R = 0.5, and R — 0.9, 
respectively. Other parameters as in Fig. [7] 
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FIG. 10: (Color online) Power spectrum of the summarized signal xs ~ xi + X2 for a fixed time delay of r = 5. Panels (a), 
(b), and (c) correspond the case of moderate (C = 0.2, Di = 0.6), weak (C = 0.1, Di = 0.6), and strong (C — 0.2, Di — 0.15) 
synchronization. The solid (black), dashed (red), and dotted (green) curves refer to a memory parameter of i? = 0, 0.5, and 
0.9, respectively. Other parameters as in Fig. [T] 



ate a phase from the time series of spikes: 

ifit) = 27r ^~^"^ +2^(^-l), 



where ti denotes the time of the z-th spike. With 
this definition the phase increases by a value of 27r 
for each spike 



iPikovskv et al 



Hauschildtl. 120051: iHauschildt et di [200I 
1996ll200l| . The phase difference A(p 



between two subsystems can be defined for general n : m 
synchronization as follows 



\(pi{t) V2(i)|, 

n 



where ipi(t) and f2it) denote the phases of the respec- 
tive subunits. In this work, we consider only 1 : 1- 

synchronizat ion. 

In Refs. [Hauschildt et al. I, I2OO6I: iRosenblum et all 
a measure for the phase synchronization is consid- 
ered: the so-called synchronization index 7. This quan- 
tity is also defined using the phase difference Aip: 

7 = ^(cosA(p(t))2 + (sinA(/3(t))2. 



It varies between (no synchronization) and 1 (perfect 
synchronization) . 

(4) If the two subsystems are already strongly synchro- 
nized, it is helpful to consider the time intervals dur- 
ing which the phase difference stays in a 27r-phase ra nge 
jLai et aLl . [2OO6I : iPark fc Lai . 120051 [Park et al\ . \^Ol\ . In 
these synchronization intervals, the subsystems exhibit 
the same number of spikes. In the case of strong synchro- 
nization, the synchronization index 7 shows only small 
changes near its maximum value, whereas the time inter- 
val of constant phase difference can vary significantly. 

(5) A measure of the amount of synchronization is given 
by the average length S of the synchronization intervals. 
The stronger sensitivity, as compared to 7, can be seen in 
Fig. I16[ which depicts the average phase synchronization 
interval S as black dots (solid curve) and the synchro- 
nization index 7 as red squares (dashed line) for varying 
time delay t. The other control parameters are chosen 
as i? = and K = 1.5. The scale is chosen such that the 
points for r = and r = 10 coincide. Note that the av- 
erage phase synchronization interval S shows a stronger 

(6) increase than the synchronization index 7. Since the av- 
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FIG. 11: (Color online) Interspike interval (Te) distribution of the summarized signal = Xi+X2 for moderate synchronization 
(C = 0.2, Di — 0.6) in dependence on the time delay r. The greyscale (color coding) denotes the probability of finding a 
certain interspike interval Te- Panels (a), (b), and (c) correspond to a memory parameter of R = 0, R — 0.5, and R — 0.9, 
respectively. The white and grey dotted lines in panel (a) at r = Te and r — Ts/2, respectively, are guides to the eye. Other 
parameters as in Fig. [7] 




FIG. 12: (Color online) Interspike interval (Te) distribution of the summarized signal xs = xi + X2 for weak synchronization 
(C = 0.1, Di = 0.6) in dependence on the time delay r. The greyscale (color coding) denotes the probability of finding a 
certain interspike interval Te. Panels (a), (b), and (c) correspond to a memory parameter oi R = Q, R = 0.5, and R = 0.9, 
respectively. Other parameters as in Fig. [T] 



erage phase synchronization interval S is more sensitive 
to the control parameters K, r, and R than the synchro- 
nization index 7, we restrict om' investigations to the 
discussion of S. 

Figure [12] displays the time evolution of the phase dif- 
ference Aip, when no control is applied to the system. 
The solid (black), dashed (red), and dash-dotted (green) 
curves correspond to the case of moderate (C = 0.2, 
Di = 0.6), weak (C = 0.1, Di = 0.6), and strong 
(C = 0.2, Di — 0.15) synchronization, respectively. The 
inset depicts an enlargement for moderate synchroniza- 
tion. From this inset, one can clearly see the plateaus in 
between phase jumps of 27r. At these jumps, only one 
subsystem shows a spike whereas the other one remains 
subthreshold. For better synchronization, the slope of 
Aif becomes flatter. In the case of strong synchroniza- 
tion, for instance, Aip shows only a few phase jumps and 
remains in a 27r-range for large time intervals. The quan- 



tity S measures the average lengths of these intervals. 

Figure flTl shows effects of extended time-delayed feed- 
back in the average phase synchronization interval S for 
varying time delay r. The feedback gain K is fixed at 
K = 1.5. Panels (a), (b), and (c) display the case of 
moderately, weakly, and strongly synchronized subsys- 
tems, respectively. In all panels, the solid (black), dashed 
(red), and dotted (green) curves refer to a memory pa- 
rameter of i? = 0, i? = 0.5, and R = 0.9, respectively. 
The insets in panels (a) and (b) are enlargements for 
small time delays and the inset in panel (c) display large 
values of 6. 

In general, time-delayed feedback enlarges the average 
phase synchronization interval 5. Especially for small 
time delays, e.g., r = 0.7 for R = 0, S becomes substan- 
tially larger. For i? = 0, a modulation of 5 can be seen 
for small delays, see insets in all panels of Fig. [T71 These 
deviations are less pronounced for increasing R. Only 
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FIG. 13: (Color online) Interspike interval (Te) distribution of the summarized signal xj: — xi + X2 for strong synchronization 
(C = 0.2, Di — 0.15) in dependence on the time delay r. The greyscale (color coding) denotes the probability of finding a 
certain interspike interval Te- Panels (a), (b), and (c) correspond to a memory parameter of i? = 0, i? = 0.5, and R = 0.9, 
respectively. Other parameters as in Fig. [7] 



(a) moderate (b) weak (c) strong 
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FIG. 14: (Color online) Interspike interval (Te) distribution of the summarized signal xe — xi+X2 in dependence on the time 
delay r for a memory parameter 7? — 0.99 for moderate (C — 0.2, Di = 0.6), weak {C — 0.1, -Di = 0.6), and strong (C — 0.2, 
Di = 0.15) synchronization in panels (a), (b), and (c), respectively. The greyscale (color coding) denotes the probability of 
finding a certain interspike interval Te. Other parameters as in Fig. [7] 



panel (c), which refers to strong synchronization, shows 
larger values of S with increasing memory parameter. 

The sensitivity of S, as discussed in Fig. [161 can also be 
seen in the case of strong synchronization. See panel (c) 
in Fig. 1171 Since the already strong synchronization is 
further enhanced by the control force, the average phase 
synchronization interval rises by several orders of magni- 
tude as shown in the inset. For perfect synchronization 
and simultaneous spiking, S would be arbitrarily large 
and merely reflect the integration time. 



7. CONCLUSION 

We have investigated the cooperative dynamics of two 
symmetrically coupled neurons under the effects of local 
extended time-delayed feedback. We have found that the 
specific choice of the control parameters, i.e., feedback 



gain K, time delay r, and memory parameter i?, alters 
the cooperativity, which we have discussed in the context 
of different measures of synchronization: ratio of average 
interspike intervals and average phase synchronization in- 
tervals. If the control force is generated including states 
further in the past, i.e., for larger memory parameter, the 
two subsystems exhibit enhanced phase synchronization, 
if the uncontrolled system is already in the strong syn- 
chronization regime. In general, both the frequency and 
phase synchronization as well as the interspike interval 
distribution become less sensitive to variations in r with 
increasing R. A small memory parameter R leads to a 
suppression of interspike intervals slightly larger than the 
time delay. 

The stability of a synchronous manifold for coupled 
ordinary differential equations can be investigated by 
means of a master stabi lity function [Chavez et aLl . [2005t 
iPecora fc Carroll Il998| . For delay differential equations 
involving stochastic input, this formalism is not yet devel- 
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oped. Future investigations should address this to obtain 
an analytical understanding of the control of synchroniza- 
tion in neural networks by time-delayed feedback. 



FIG. 15: (Color online) Phase difference in units of 2n. The 
solid (black), dashed (red), and dash-dotted (green) curves 
correspond to the case of moderate (C — 0.2, Di = 0.6), 
weak (C = 0.1, Di = 0.6), and strong (C = 0.2, Di = 0.15) 
synchronization, respectively. No control is applied to the 
system. The inset shows an enlargement for moderate syn- 
chronization. Other parameters as in Fig. (2] 
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FIG. 16: (Color online) Average phase synchronization inter- 
val S (black dots, solid curve) and synchronization index 7 
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